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VORTEX  DYNAMICS  OF  COHERENT  AND  CHAOTIC  STRUCTURES 
(Including  Algorithms  for  Computer  Simulations 
and  Diagnosis) 

Proposal  No.  21371-MA 

During  this  three  year  period  we  made  substantial  contributions  to  the  understanding  of  fundamental 
processes  in  two-dimensional  in  viscid  and  nearly  dissipationless  vortex  dynamics  through  an  approach  we  call 
complementary  modeling.  We  have  identified  and  studied  physical -space  processes  like  merger,  binding, 
axisymmetrization-and-gradient-intensification  of  near-isolated  regions  of  vorticity.  We  have  employed  a  variety 
of  algorithms  and  codes  including  pseudo-spectral,  contour  dynamical,  and  our  recently  introduced  moment 
model.  Through  the  last  application,  we  have  solved  analytically  an  asymptotical  representation  of  the  sym¬ 
metric  merger  problem  and  provided  causes  and  conditions  for  merger  that  agree  with  pseudospectral  and  con¬ 
tour  dynamical  simulations.  We  have  also  delineated  the  essential  parameters  which  control  the  final  state  of 
two  nonidentical  nearby  vortices  (asymmetric  merger) 

We  have  established  an  Interactive  Graphics  and  Numerical  Diagnostics  Laboratory  to  examine  data  from 
large-scale  computer  simulations.  This  uses  SUN  Microsystems  and  Hitachi  hardware  (SUN  models  2/120). 
This  allows  us  to  diagnose  or  "project"  large  quantities  of  field  data  like  vorticity  or  velocity,  into  lower  dimen¬ 
sional  forms.  The  system  will  accelerate  the  creative  process  and  will  enhance  the  productivity  of  scientists  and 
engineers  in  computational  science. 

As  a  result  of  interactions  with  government  laboratory  personnel  at  the  ARO  sponsored  Blade-Vortex 
interaction  workshop  at  NASA  Ames,  we  began  a  consideration  of  compressible  and  3D  vortex  dynamics.  Work 
with  F.  Lund  presented  a  new  formalism  to  treat  perturbation  terms  normally  omitted  from  "aeroacoustic”  theory 
which  give  the  effect  of  sound  on  vortex  motions.  This  was  made  explicit  in  3D  for  the  local  induction  approxi¬ 
mation.  M.  V.  Melander  took  the  lead  in  examining  the  3D  incompressible  evolution  of  two  offset-orthogonal 
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Gaussian  vortex  distributions,  namely,  the  "reconnection"  problem,  a  quintessential  process  in  3D  turbulence  and 
one  which  is  probably  the  chief  cause  of  intermittency.  His  low  dissipation  643,  963,  and  1283  dealiased 
spectral  simulations  iof  the  Navier-Stokes  equations  show  an  inviscid  entanglement  on  a  convective  time  scale 
which  yields  an  apparent  visual  reconnection  when  the  vortieity  distributions  are  "coarse  grained."  The  essential 
interactions  of  "hairpin" -like  structures  pulled  orthogonally  from  the  original  cores  and  "spikes"  ejected  from  the 
close  binding  of  oppositely  directed  vortex  tubes  are  features  that  dominate  die  inviscid  entanglement  process.. 

Toward  the  end  we  initiated  a  collaboration  with  K.-H.  Winkler  of  Los  Alamos  and  P.  L.  Woodward  of  die 
University  of  Minnesota.  We  have  done  simulations  at  Los  Alamos  and  used  their  ultra-speed  graphics  facility. 
Using  Woodward’s  PPM  code,  we  have  simulated  the  interaction  of  an  M  -  1.25,  2.5  and  3.0  shock  with  heavy 
and  light  bubbles,  in  2D  axisymmetry  (z^j,  and  overlapped  the  parameter  region  of  the  Haas-Sturtevant  (Cal¬ 
tech)  experiments.  The  break-up  of  bubbles,  primarily  by  vortex  deposition  at  the  interface  is  under  investiga¬ 
tion.  A  paper  on  the  subject  was  submitted  to  the  Tokyo  meeting  on  Fundamentals  of  Vortex  Motion  and  a 
prelimn  ary  account  was  published  in  "A  Numerical  Laboratory,"  Phys.  Today,  October  1987.  This  work  has 
potential  relevance  for  the  multiple  nuclear  burst  in  a  dust  (inhomogeneity)  filled  environment 

We  believe  the  last  three  areas  represent  virgin  territory  in  computational  and  mathematical  fluid  dynamics 
and  are  relevant  to  many  DOD  applications.  We  hope  that  funds  will  be  available  to  us  soon  to  pursue  these 


R 


Kb 


f«n V*  V  V *  -  .  ,♦  ;  V  -  V V  »Vr V* ' 


REFERENCES 


1.  Zabusky,  N.  J.,  Computational  synergetics  and  the  exploration  of  nonlinear  science.  Lectures  in  Mathemati¬ 
cal  Physics  10  (1985)  143-147. 

2.  Zabusky,  N.  J.t  Visualizing  Mathematics.  Evolution  of  vortical  flows,  Physica  18D  (1986)  15-25. 

3.  Melander,  M.  V.,  Zabusky,  N.  J.  and  Styczek,  A  moment  model  for  vortex  interactions  of  the  2D  Euler 
equations,  J.  Fluid  Mech.,  Vol.  167  (1986),  95-115. 

4.  Melander,  M.  V.,  Overman,  E.  A.  and  Zabusky,  N.  J.,  Computational  vortex  dynamics  in  two  and  three 
dimensions,  Appld.  Numerical  Mathematics  3  (1,2)  (1987),  59-80. 

5.  Melander,  M.  V.,  McWilliams,  J.  C.  and  Zabusky,  N.  J.,  Axisymmetrization  and  vorticity-gradient 
intensification  of  an  isolated  two-dimension  vortex  through  filamentation,  J.  Fluid  Mech.  178  (1987)  137- 
159. 

6.  Lund,  F.  and  Zabusky,  N.  J.,  A  new  approach  to  the  effect  of  sound  on  vortex  dynamics,  Phys.  Fluids  30 
(1987)  2306,  2307. 

7.  Q.  Zou,  Overman,  E.  A.,  II,  Wu,  H.-M.  and  Zabusky,  N.  J.,  Contour  dynamics  for  the  Euler  equations. 
Curvature  controlled  initial  node  placement  and  accuracy,  J.  Comput.  Physics,  in  press. 

8.  Melander,  M.  V.,  Zabusky,  N.  J.  and  McWilliams,  J.  C.,  Symmetric  vortex  merger  in  two-dimensions: 
Causes  and  conditions,  J.  Fluid  Mech.,  in  press. 

9.  Seren,  C,  Melander,  M.  V.  and  Zabusky,  N.  J.,  Phys.  Fluids  30  (1987)  2604.  Color  Montage,  winner  at  the 
Picture  Gallery  of  Fluid  Motions,  APS  DFD  meeting  in  Columbus,  Ohio,  1 1/86. 

10.  Melander,  M.  V.  and  Zabusky,  N.  J.,  Entanglement  produces  an  apparent  reconnection  of  3D-vortex  tubes, 
Phys.  Rev.  Letters,  submitted  12/87. 

11  Winkler,  K.-H.,  Chalmers,  J.  W„  Hodson,  S.  W,  Woodward,  P.  R.  and  Zabusky,  N.  J.,  A  Numerical 
Laboratory,  Physics  Today  (1987),  28-37. 


i 


1 

) 


J 


